Breakdown of thermalization in finite one-dimensional systems 
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We use quantum quenches to study the dynamics and thermalization of hardcore bosons in finite one- 
dimensional lattices. We perform exact diagonalizations and find that, far away from integrability, few-body 
observables thermalize. We then study the breakdown of thermalization as one approaches an integrable point. 
This is found to be a smooth process in which the predictions of standard statistical mechanics continuously 
worsen as the system moves toward integrability. We establish a direct connection between the presence or 
absence of thermalization and the validity or failure of the eigenstate thermalization hypothesis, respectively. 
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Little more than fifty years ago, Fermi, Pasta, and Ulam 
(FPU) set up a numerical experiment to prove the ergodic 
hypothesis for a one-dimensional (ID) lattice of harmonic os- 
cillators once nonlinear couplings were added. Much to their 
surprise, the system exhibited long-time periodic dynamics 
with no signals of ergodicity. This behavior could not be 
explained in terms of Poincare recurrences and motivated in- 
tense research |2], which ultimately gave rise to the modern 
chaos theory. It led to the discovery of solitons (stable soli- 
tary waves) in nonlinear systems and to the understanding of 
thermalization in terms of dynamical chaos [2]. In the lat- 
ter scenario, there is a threshold below which the interactions 
breaking integrability are ineffective in producing chaotic be- 
havior and the system cannot be described by standard statisti- 
cal mechanics ||3|]. The FPU numerical calculations happened 
to be below that threshold fl. 

More recently, experiments with ultracold gases in ID ge- 
ometries have challenged our understanding of the quantum 
domain fl. After bringing a nearly isolated system out of 
equilibrium, no signals of relaxation to the expected thermal 
equilibrium distribution were observed. Some insight can be 
gained in the framework of integrable quantum systems |@], 
but then it remains the question of why thermalization did 
not occur even when the system was supposed to be away 
from integrability. In the latter regime, thermalization is ex- 
pected to occur f?', 's'l. This new experimental result (E] has 
opened many questions such as: Will thermalization occur if 
one waits longer? Is there a threshold after which thermal- 
ization will occur? In this work we address some of these 
questions using numerical experiments. 

In the limit in which the quantum system is integrable, it has 
been shown numerically 161] that observables such as the ones 
measured experimentally relax to an equilibrium distribution 
different from the thermal one. That a novel distribution is 
generated is the result of the conserved quantities that render 
the system integrable, and can be characterized by a general- 
ization of the Gibbs ensemble (GGE) [6]. Several works since 
then have addressed the relevance and limitations of the GGE 
to various integrable systems and classes of observables 
Much less is known away from integrability where fewer ana- 
lytical tools are available and numerical computations become 
more demanding. Early works in ID have provided mixed re- 



sults; thermalization was observed in some regimes and not 
in others [10]. In two dimensions, thermalization was unam- 
biguously shown to occur iQ] and could be understood on the 
basis of the eigenstate thermalization hypothesis [f\. Recent 
works have also pointed out a possible intermediate quasi- 
steady regime that could occur before thermalization in a class 
of fermionic systems fll']. Here we study how breaking inte- 
grability affects the thermalization of correlated bosons in a 
ID lattice after a quantum quench. 

We consider impenetrable bosons in a periodic ID lattice 
with nearest-neighbor hopping t and repulsive interaction V, 
and next-nearest-neighbor hopping t' and repulsive interac- 
tion V'. The Hamiltonian reads 11211 



-t' (blb,+2 + ii.c 



(1) 



When t' = V' = this model is integrable. In order to 
understand how the proximity to the integrable point affects 
equilibration, we prepare an initial state that is an eigenstate of 
a system with t — tini, V = Vini, t', V and then quench the 
nearest-neighbor parameters to t = tfin, V = Vfin without 
changing t', V, i.e., we only change Um,Vini tf^n, Vf^n■ 
The same quench is then repeated for different values of t', V' 
as one approaches t' = V' — 0. We have performed the exact 
time evolution of up to eight impenetrable bosons in lattices 
with up to 24 sites. Taking advantage of translational invari- 
ance, this required the full diagonalization of blocks in the 
Hamiltonian that contained up to 30,667 states. 

Does integrability, or its absence, affect the relaxation dy- 
namics of experimentally relevant observables? To answer 
that question, we examine two of those observables: the mo- 
mentum distribution function n{k) and the structure factor for 
the density-density correlations N{k) lfl2ll . Since the initial 
state wavefunction can be expanded in the eigenstate basis of 
the final Hamiltonian H as \ipini} — J2a CqI^q), one finds 
that, if the spectrum is nondegenerate and incommensurate, 
the infinite time average of an observable O can be written as 



{O) — Odiag — ^ ^ \Ca\^Oaa, 



(2) 



2 



where Oaa are the matrix elements of O in the basis of the 
final Hamiltonian. This exact result can be thought as the pre- 
diction of a "diagonal ensemble", where |Cq,|^ is the weight 
of each state of this ensemble |8]. We then study the normal- 
ized area between our observables, during the time evolution, 
and their infinite time average, i.e., at each time we compute 



5nk = 



T,k Hk) - nd,ag{k)\ 



(3) 



and similarly for 5Nk- If, up to small fluctuations, n{k) 
and N{k) relax to a constant distribution, it must be the 
one predicted by the infinite time average in Eq. dU, i.e., 

5n{k), 5N{k) 0. 
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FIG. 1: Quantum quench from tini = 0.5, Vini = 2.0tof/i„ = 1.0, 
V/i„ = 1.0, with f-„i = t)i^ = t' and V^i = V/i„ = V' in a 
system with eight hardcore bosons {N\, = 8) in 24 lattice sites (L = 
24). The initial state was chosen within the eigenstates of the initial 
Hamiltonian with total momentum A: = in such a way that after the 
quench the system has an effective temperature T = 2.0 in all cases. 
Given the energy of the initial state in the final Hamiltonian E = 
{il)ini\Hfiri\il}i„i) the effective temperature is computed as described 
in Ref. (131]. (a) Initial and diagonal ensemble results for n{k) when 
t' ^ 0.24. (b)-(e) Time evolution of Suk for t' = V = 0.24, 

0.12, 0.03, and 0.0 respectively, (f) Initial and diagonal ensemble 
results for N{k) when t' = V' — 0.24. (g)-(j) Time evolution of 
SNk for t' ^V' = 0.24, 0.12, 0.03, and 0.0 respectively 

In Fig. [1] we show results for 6nk and SNk as a function 
of time T for four different quenches as one approaches the 
integrable point. In Figs.llla) and[Tlf), we compare the initial 
n{k) and N{k) with the predictions of Eq. away from in- 
tegrability. Interestingly, the time evolution of 6nk and SNk 
in Figs.[nb)-[TIe) and Figs.[ng)-[Ilj), respectively, can be seen 
to be very similar for all values of t' and V'. Hence, the dy- 
namics are barely affected by the closeness to the integrable 
point Iil2i] . In all cases, we find that there is a fast relaxation 
of n{k) and N{k) towards the diagonal ensemble prediction 
(in a time scale t ~ 1/tfin Ill2ll) and that later they fluctu- 
ate within one to two percent of that prediction. The latter 
differences decrease as one increases the system size and the 
energy, or the effective temperature jisll of the system fl^. 
From these results, we infer that, for very large systems sizes 



and finite effective temperatures, n{k) and N{k) should in 
general relax to exactly the predictions of Eq. (|2|i even if the 
system is very close or at integrability. 

Once it is known that the diagonal ensemble provides a 
very good description of observables after their relaxation dy- 
namics, a question that remains to be answered is how good 
are standard statistical ensembles in reproducing the n{k) and 
N{k) distributions predicted by the diagonal ensemble as one 
approaches the integrable point? Based on the experimental 
results (jstl one would expect conventional statistical mechan- 
ics to fail everywhere in ID or at least over a finite window in 
the vicinity of the integrable point. We find the latter scenario 
to be the one realized in our finite ID lattices. 

In the main panel of Fig. |2la), we compare the diagonal 
ensemble results with the predictions of the microcanonical 
ensemble for our two observables of interest. We find the 
differences between them to be below one percent and de- 
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FIG. 2: Comparison between the predictions of different statisti- 
cal ensembles for n{k) and N{k) after relaxation in a system with 
i = 24 and Nb = 8. Results are shown for T = 2.0 and T = 3.0. 
(a) (main panel) Relative difference between the predictions of the 
diagonal ensemble and the microcanonical ensemble, (a) (inset) Rel- 
ative difference between the predictions of two diagonal ensembles 
generated by different initial states selected from the eigenstates of 
a Hamiltonian with tini = 0.5, Vini ~ 2.0 and a Hamiltonian with 
tini = 2.0, Vini = 0.5, and L = 24, Nt = 8. The final Hamil- 
tonian (tfin = 1.0, Vfin ~ 1.0) and the effective temperature are 
identical for both initial states. Relative differences are computed 
following Eq. (b) Average relative deviation of eigenstate expec- 
tation values with respect to the microcanonical prediction (see text) 
as a function of t' ,V' , for the same effective temperatures as in (a). 
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creasing with system size lll2|] when the system is far from 
integrabiUty (t' = V' > 0.1). Thus, one can say that ther- 
maHzation takes place in this case. As one approaches the 
integrable point, on the other hand, the differences between 
the diagonal and microcanonical ensembles increase, signal- 
ing a breakdown of thermalization in ID. This breakdown is 
accompanied by a dependence of the properties of the system 
after relaxation on the initial state, as can be seen in the inset 
in Fig. 12 a). There, we compare the predictions of the diago- 
nal ensemble for two different initial states that have the same 
effective temperature. Those predictions are almost indepen- 
dent of the initial state for t' — V' > 0.1, but strongly depen- 
dent on it closer to integrability. In Fig.|2a), it is remarkable 
to find that both n{k) and N{k) exhibit very similar quantita- 
tive behavior away from integrability, which is an indication 
of the generality of our results for few-body observables. 

A clarification is in order at this point. It is usually assumed 
that, for extensive quantities, the predictions of conventional 
statistical mechanical ensembles such as the microcanonical 
and canonical ones are identical, provided the systems are 
large enough. We should stress that this is not the case in our 
small systems, and may not be the case for many experiments 
that are performed in similar setups. In Fig.[3]one can see that 
depending on the physical observable under consideration the 
differences between the predictions of those ensembles can 
be quite large [in particular for n{k)]. Interestingly, they are 
almost not affected by the proximity to integrability, and, as 
expected, can be seen to reduce with increasing system size. 
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FIG. 3: Relative difference between the predictions of the micro- 
canonical and canonical ensembles for our two observables of inter- 
est in systems with L = 21, A^j, = 7 and L = 24, iV^ = 8. Relative 
differences are computed following Eq. Jsj. 

One may wonder why conventional statistical mechanical 
ensembles can predict the exact values of few-body observ- 
ables after relaxation at all, since these seem to be a priori 
fully dependent on the initial conditions through the values of 
Cq p [see Eq. ©I. In the main panel in Fig. SJb), we show 
the distribution of |Cq p for one of the quenches that exhibited 
thermalization in Fig.|2fa). FigureHfb) shows that the |CqP 
distribution is clearly different from the microcanonical one, 
in which all states within an energy window are taken with the 
same weight. One can conclude then that something else is at 
play here. 



A resolution to this puzzle was advanced by Deutsch and 
Srednicki in terms of the eigenstate thermalization hypothe- 
sis (ETH) (01 • ETH states that the fluctuation of eigenstate 
expectation values of generic few-body observables [Oqq in 
Eq. ©I is small between eigenstates that are close in energy, 
which means that the microcanonical average is identical to 
the prediction of each eigenstate, or what is the same that the 
eigenstates aheady exhibit a thermal value of the observable. 
If this holds, thermalization in an isolated quantum system 
will follow for any distribution of |CqP, as long as it is nar- 
row enough in energy. This scenario was shown to be valid in 
isolated two-dimensional systems in Ref. |8]. 
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FIG. 4: (a) n(fc — 0) as a function of the energy for all the eigenstates 
of the Hamiltonian (including all momentum sectors), (main panel) 
t=V ^landt' = V' ^ 0.24. (inset) f = V = 1 and t' = V = 
0.03. The system has 24 lattice sites and 8 bosons for which the total 
Hilbert space consists of 735,471 states, (b) Distribution of |Cc, ^ for 
the quench from tini — 0.5, Vini = 2.0 to tfin = 1.0, V/i„ — 1.0 
for an effective temperature T — 2.0. (main panel) t = V = 1 and 
t' = V' ^ 0.24, where E = {ip^n^\Hf;„\ip^n^) = -6.85. (inset) 
t = V = lundt' = V' = 0.03, where E = -5.95. 



In ID systems, we find the onset of thermalization to be 
directly related to the validity of the eigenstate thermalization 
hypothesis. In the main panel in Fig.Ua), we depict n{k = 0) 
[similar results were obtained for n{k ^ 0) and iV(fc)] in each 
eigenstate of the Hamiltonian as a function of the energy of the 
eigenstate, when the system is far from integrability. Figure 
HJa) shows that after a region of low energies where the eigen- 
state expectation values exhibit large fluctuations follows an- 
other region where fluctuations are small and ETH holds. The 
inset in Fig. HJa) shows that for a system close to the inte- 
grable point, in which thermalization is absent [Fig.|2fa)], the 
eigenstate to eigenstate fluctuations of n(fc = 0) are very large 
over the entire spectrum and ETH does not hold. 

In order to be more quantitative, and to understand how 
ETH breaks down as one approaches the integrable point, we 
have computed the average relative deviation of the eigenstate 
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expectation values with respect to the microcanonical predic- 
tion (A™'^) for n{k = 0) and for N{k = tt) [12]. They 
are shown in Fig. ^a) and|5fb), respectively, as a function of 
the effective temperature T corresponding to an energy E of 
the microcanonical ensemble [13]. Temperature in this plot is 
only used as an auxiliary tool for assessing how far from the 
ground state these systems are. 

In Fig.|5]one can see that A™"=n(fc = 0) and A""=iV(fc = 
tt) exhibit very similar behavior Below T ~ 1.5 fluctuations 
are large and due to the relatively small number of states in 
some energy windows the use of the microcanonical ensem- 
ble may not be well justified. For T > 1.5, the fluctuations 
saturate with increasing temperature, and A™*'^n(fc = 0) and 
A'"*'^A^(fc — tt) continuously increase, at any given temper- 
ature, as one approaches the integrable point. The latter is 
better seen in Fig. |2b) for the two specific effective temper- 
atures studied in Fig. EJa). A comparison between Fig. |2ja) 
and Fig. |2b), allows us to stablish a direct connection be- 
tween the breakdown of thermalization in ID and the increase 
of the eigenstate to eigenstate fluctuations in the Hamiltonian, 
or what is the same, the failure of ETH. We have also stud- 
ied how these results are modified when changing the system 
size. We find that for our lattice sizes A™"^ still reduces as 



one increases the system size 111211 . The relative reduction is 
more pronounced away from integrability. 
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FIG. 5: Average relative deviation of eigenstate expectation values 
with respect to the microcanonical prediction as a function of the 
effective temperature of the microcanonical ensemble. Results are 
presented for: (a) n{k = 0) and (b) N{k — tt). These systems have 
24 lattice sites, 8 bosons, and t — V — 1.0. An effective temperature 
T = 2 in these plots corresponds to an energy Emic ~ —5.95 for 
t' = V' ^ 0.03 and Er^ic = -6.85 for t' = V' ^ 0.24 and, T = 
10 corresponds to an energy Emic ~ —0.81 for t' — V' — 0.03 and 
Emic = —0.82 for t' = V' = 0.24 (see the corresponding parts of 
the spectrum in Fig.|4l(. 

Overall, our study shows that in finite ID systems there is a 
regime close to integrability where thermalization fails to oc- 
cur, and that this failure originates in the breakdown of ETH. 
We also find that ETH holds further departing from integrabil- 
ity explainingwhy thermaUzation can occur in ID. Since both 
experiments [5] and computations have been performed with 
relatively small systems, an important question that remains 
open is whether for sufficiently large system sizes thermal- 
ization will occur arbitrarily close to the integrable point or 



whether there will be a finite critical value starting from which 
integrability breaking terms will produce thermalization. For 
our finite systems all we see is a smooth breakdown of ther- 
malization as one approaches integrability |14]. The answer 
to this question may be well suited for experimental analysis 
as it may be easier to study the scaling with system size in 
experiments than within our numerical computations, which 
scale exponentially with the system size. The addition of a 
lattice along the ID tubes in the experiments iflill may help in 
addressing these issues. 
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town University and by the US Office of Naval Research 
Award No. NOOO 1409 10966. We are grateful to Maxim Ol- 
shanii for discussions and useful comments on the manuscript. 
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Hamiltonian and Numerical Calculations 

In a system of units where h ~ 1, the Hamiltonian for 
our model of impenetrable bosons (hardcore bosons) in a one- 
dimensional lattice with periodic boundary conditions reads 



-t'iUh 



1+2 



H.c, 



H.c 



v 



V 



ni — - 
2 



(4) 



where the hardcore boson creation and annihilation operators 
at site i are denoted by and Sj, respectively, and the local 
density operator by hj = blb^. In different sites the creation 
and annihilation operators commute as usual for bosons: 

;;t if 



Kb]] 



0, for i 7^ j. 



(5) 



However, on the same site the hardcore bosons operators sat- 
isfy anticommutation relations typical of fermions: 



1, 



Lt2 



0. 



(6) 



These constraints avoid double or higher occupancy of the lat- 
tice sites. In Eq. (|4|i, the nearest and next-nearest-neighbor 
hopping parameters are denoted by t and t', the nearest and 
next-nearest-neighbor interactions are denoted by V and V', 
which in our study are always repulsive, i.e., V, V > 0, and 
L is the number of lattice sites. 

For t' = V ^ 0, the above Hamiltonian is integrable. It 
also exhibits a superfluid to insulator transition at half-filling 
as V is increased [{V/t)c — 2]. Finite values of t' and 
V' break integrability and generate a plethora of competing 
phases. All our results in this work are obtained for quenches 
within the superfluid phase. 

For our exact study of the nonequilibrium dynamics, we 
perform a full diagonalization of the many-body Hamiltonian 
dUi taking advantage of the translational symmetry of the lat- 
tice. Our initial state is always selected from the zero total 
momentum sector. The largest fc = sector in our study cor- 
responds to eight bosons in 24 lattices sites, where matrices of 
dimension D = 30, 667 are diagonalized. Since our quenches 
do not break translational symmetry only states with fc = 
are required for the time evolution, and we use all of them 



|V(t)) = e-*^n^.m _ 

a 

where lipir)) is the time-evolving state 
state, l^f") are the eigenstates of the Hamiltonian with zero 
total momentum and energy Ea, and C, 



■(/'ini) is the initial 
nian 



The observables of interest in this work are the momentum 
distribution function of the bosons 



which is generally measured in ultracold gases experiments, 
and the density-density correlation structure factor 



N{k) = 



IE 



Since we always work at a fixed number of particles Ni,, the 
expectation value of N{k = 0) is trivially {N{k — 0)) = 
/L. We subtract this trivial constant in all our plots. 
Studying a one-body and a two-body observable allows us 
to assess the generality of our results for generic few-body 
observables in generic one-dimensional systems. 



Relaxation Dynamics after a Quantum Quench 

In Fig. 1 (main text), we have shown results for the re- 
laxation dynamics of a system with 8 bosons (TVf, = 8) in 
24 lattice sites (L = 24). The initial state was selected so 
that the time evolving system has an effective temperature 
T ~ 2.0. [Given the energy of the time evolving state in the 
final Hamiltonian, which is conserved E — {ipini\H fin\ipini) , 
the effective temperature T is defined by the expression E — 

Z^^Tr |ff/i„ exp(-i?/i„/fcsT)|, where Hjin is the final 

Hamiltonian, Z = Tt ^exp{—Hfin/kBT)^, and fc^ the 
Boltzmann constant. The trace runs over the full spectrum, 
including all momentum sectors.] In Fig.|6] we present results 
for a quench with exactly the same Hamiltonian parameters 
as in Fig. 1 (main text) but taking a different initial state so 
that the energy of the system is higher, corresponding to an 
effective temperature T = 3.0. (Notice that for all our results 
tfin — 1-0 sets the energy scale.) 

Overall, we find that the time scale over which relaxation 
occurs is of the order t ^ 1/t fin in all cases studied. The 
short time scale behavior can be better seen in the insets in 
Fig lSg)Si)- In addition, the results in Fig.|6] for an effective 
T = 3.0, show that the values of Suk and 6Nk after relaxation 
are slightly smaller and fluctuate less than the ones obtained 
for an effective temperature T = 2.0 [Fig. 1 (main text)]. This 
implies that increasing the energy of the final state improves 
the convergence towards the infinite time average prediction. 
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This difference can be understood in terms of a higher num- 
ber of eigenstates of the Hamihonian that participate dur- 
ing the nonequilibrium dynamics as the energy of the time- 
evolving state is increased. Since more states are involved, 
dephasing is more effective and time fluctuations after relax- 



ation are smaller in these finite systems. The higher partic- 
ipation can be understood as the density of states increases 
with increasing energy, so that in general there are more states 
available close to any given energy the higher the latter is. 
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FIG. 6: Quantum quench from Uni — 0.5, Vini — 2.0 to tfin — 1.0, V/i„ = 1.0, with t'^^i — t^^, 
with eight hardcore bosons (A^;, = 8) in 24 lattice sites (L = 24). The initial state was chosen within the eigenstates of the initial Hamiltonian 
with total momentum fc = in such a way that after the quench the system has an effective temperature T — 3.0 in all cases (see text), (a) 
Initial and diagonal ensemble results for n{k) and t' — V' = 0.24. (b)-(e) Time evolution of Suk for t' = V' = 0.24, 0.12, 0.03, and 0.0 
respectively, (f) Initial and diagonal ensemble results for N{k) and t' = V' — 0.24. (g)-(j) Time evolution of SNk for t' = V' = 0.24, 0.12, 
0.03, and 0.0 respectively, (g)-(i) (insets) Magnification of the short-time scales for the time evolution depicted in the main panels. 



In Fig. |7] we depict the number of states that have a value 
of I Cq I ^ greater than 10^^, 10^'', 0.1, when the effective 
temperature is 2.0 [Fig.|7ta)] and 3.0 [Fig.EJb)]. One can see 
there that the number of states with |CqP > lO"'', |CaP > 
lO^'', and |Cap > 10"^ j^rger for T = 3.0 flian for T = 
2.0 (notice that this is a log-log plot). Since = 1' 

this means that at lower temperatures there are more states 
with larger values of | | ^, but since those are very few in any 
case, dephasing tends to be less effective as the temperature 
decreases. 

The same argument that allows us to understand the effect 
of the effective temperature helps in understanding why at in- 
tegrability, or very close to it, Sn^ and SN^ exhibit larger time 
fluctuations and larger average values. The effect is small 
(~ 1%) in Fig. 1 (main text) and Fig.|6] but can be also under- 
stood by analyzing Fig.|7] There, one can see that the number 
of states with the largest values of |Cq,P is always the largest 
at integrability and continuously evolves as one moves away 
from integrability. In Fig. |7lb) such evolution is seen to sat- 
urate for t' = V' > 0.12. The effect of integrability on the 
number of states with large |Cq |^ is related to the presence of 
additional conserved quantities, which given the initial condi- 
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FIG. 7: Number of states with |Ca |^ greater than the value presented 
in the x axis, for an effective temperature (a) T — 2.0 and (b) T = 
3.0, and for exactly same quenches studied in Fig. 1 (main text) and 
Fig.m respectively. 

tions, restrict the number of states close to the energy of the 
system that can have a significant overlap with the initial state. 

We should add that the fact that in all cases time fluctuations 
are found to be relatively small, even though these systems 
are not only finite but small, can be related to the smallness of 
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the off-diagonal elements of the physical observables under 
consideration and has been discussed before yj,|2i]. 

Finally, the effect of changing the system size can be stud- 
ied by considering a smaller system with exactly the same 
density and effective temperature as the one in Fig.|6l Results 
for a system with 7 hardcore bosons in 21 lattice sites are de- 
picted in Fig.[8] They exhibit a qualitatively similar behavior 
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to the one seen in Fig. |6] but after relaxation 6nk and SNk 
show larger time fluctuations and the values around which 
they fluctuate are also larger than for the system with A^;, = 8 
and L = 24. This supports the expectation that for sufficiently 
large system sizes time fluctuations after relaxation will be 
arbitrarily small and relaxation will occur towards the exact 
predictions of the diagonal ensemble. 
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FIG. 8: Quantum quench from tini — 0.5, Vi„i — 2.0 to tfin — 1.0, Vf in = 1.0, with t'^^i — t'^^ = t' and Vl^i 



V' in a system 



with Nb = 7 and L — 21. The initial state was chosen within the eigenstates of the initial Hamiltonian with total momentum fc = in such 
a way that after the quench the system has an effective temperature T = 3.0 in all cases (see text), as was done for Fig. [6] (a) Initial and 
diagonal ensemble results for n{k) and t' — V' = 0.24. (b)-(e) Time evolution of Suk for t' = V' — 0.24, 0.12, 0.03, and 0.0 respectively, 
(f) Initial and diagonal ensemble results for N{k) and t' ^ V' = 0.24. (g)-(j) Time evolution of SNk for t' = V' ^ 0.24, 0.12, 0.03, and 0.0 
respectively. 



The Microcanonical Ensemble 



The computations of the microcanonical ensemble predic- 
tions are done averaging over all eigenstates (from all momen- 
tum sectors) that he within a window [E — AE, E + Ai?], 
where E = {tpini\H fin\tpini) , and we have taken AE — 0.1 
in all cases. We have checked that our results are robust in 
the neighborhood of this value of AE. In Fig. |9] we show 
(a) the momentum distribution function and (b) the density- 
density structure factor for three different values of AE. The 
results for each observable are indistinguishable from each 
other Further confirmation of this can be obtained looking 
at the inset in Fig. lUa), where we depict the behavior of 
n{k = 0) and N{k = tt) vs AE for 0.01 < A£; < 0.2. 
Both quantities can be seen to be almost independent of AE 
on that window. Similar results have been obtained for all 
other values of t', V' considered in this work. 
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FIG. 9: (a) (main panel) Microcanonical momentum distribution 
function and (b) density-density structure factor for three different 
values of AE. (b) (inset) n{k = 0) and N(k = tt) vs AE. The 
system has Nt = 8, L — 24, and a total energy E — —3.87, which 
corresponds to an effective temperature T — 3.0 for t' — V' = 0.03 
andt = V = 1.0. 
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Different Ensembles 

In Fig. 2(a) (main text), we compared the predictions of dif- 
ferent ensembles for our two observables of interest at differ- 
ent effective temperatures, but the same system size A^^ = 8 
and L — 24. In order to gain an understanding of finite size 
effects, in Fig.[TO]we present results for identical calculations 
in smaller systems with = 7 and L = 21, which have the 
same density as the ones in Fig. 2 (main text). Overall, we 
find a similar qualitative behavior in these smaller systems. 
However, as expected, in the smaller lattices the relative dif- 
ferences are always larger than in the larger systems. This 
is particularly evident away from integrability, where in Fig. 
2(a) (main text) all differences are below one percent while in 
Fig.[TO]they are two to three times larger than that. 




t;v' 

FIG. 10: Comparison between the predictions of different statistical 
ensembles for n{k) and N{k) after relaxation. Results are shown 
for two different effective temperatures, (a) (main panel) Relative 
difference between the predictions of the diagonal ensemble and the 
microcanonical ensemble (L — 21, Nb = 7). (a) (inset) Relative dif- 
ference between the predictions of two diagonal ensembles generated 
by different initial states, selected from the eigenstates of a Hamilto- 
nian with tini = 0.5, Vini = 2.0 and a Hamiltonian with tini = 2.0, 
Vini = 0.5. The final Hamiltonian (tfin = 1.0, V/ m = 1.0) and the 
effective temperature are identical for both initial states. 



Eigenstate Thermalization Hypotesis 

In Fig. [TlT a). we show the eigenstate expectation values of 
n{k = 0) as a function of the energy of the eigenstate for 
a system with L — 21 and Ni, — 7. The values of t, V, t' 
and V' in the main panel and the inset correspond to the exact 
same values presented in Fig. 4(a) (main text) for L = 24 and 
Nh — 8. The behavior of n{k = 0) is very similar in both 
systems in which the lattice size and the total filling are the 
only difference, as we have kept the density constant. These 
results exemplify the robustness of ETH or its failure for our 
finite systems. Notice that in Fig. 4 (main text) and Fig. (TT] 
we are showing results for the eigenstate expectation values 
of n{k = 0) in the full spectrum (including all momentum 
sectors). The validity of ETH in the full expectrum explains 
why our microcanonical calculations, in which we included 



the contributions from all momentum sectors, provide a cor- 
rect prediction for the outcome of the diagonal ensemble even 
though in the latter ensemble only states with zero total mo- 
mentum have a finite weight. 
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FIG. 11: (a) n(k = 0) as a function of the energy for all the eigen- 
states of the Hamiltonian (from all momentum sectors), (main panel) 
t=V ^l&ndt' = V' ^ 0.24. (inset) f = 1/ = 1 and = = 
0.03. The system has 21 lattice sites and 7 bosons for which the total 
Hilbert consists of 116,280 states, (b) Distribution of |Cap for the 
quench from tini = 0.5, Vini = 2.0 to tfin = 1-0, V/i„ = 1.0 
for an effective temperature T = 3.0. (main panel) t — V — 1 and 
t' = V' = 0.24, where E = {iP,m\Hf,n\tp^m) = -3.85. (inset) 
t = V = landt' = V' = 0.03, where E = {ipini\H finli^^m) = 
-3.44. 

In Fig. fTTl b). we show the distribution of | Cq p for the same 
Hamiltonian parameters as in Fig. 4(b) (main text), but once 
again for a smaller system with L = 21 and Nf, = 7, and also 
for a higher effective temperature of the time evolving state. 
In Fig. 4(b) (main text), we showed results for T = 2.0 and in 
Fig. fTTT b). we are showing results for T — 3.0. Qualitatively, 
one can see that there is not much difference between those 
two cases. The distribution of |Cq|^ is in stark contrast with 
the weight given to each eigenstate expectation value in the 
microcanonical and canonical distributions. |Cq,P vs Ea is 
a distribution that strongly depends on the initial conditions 
of the system, and |Ca|^ exhibits strong fluctuations between 
contiguous eigenstates. The only generic feature one can see 
by comparing Fig. 4(b) (main text) with Fig. [TTT b) is that in 
both cases the distribution of | Cq, p is peaked around the mean 
energy of the time evolving state. 

Deviations from ETH 

To have a quantitative understanding of how ETH breaks 
down as one approaches integrability, we computed the aver- 
age relative deviation of the eigenstate expectation values with 
respect to the microcanonical prediction, A'"*'^n(fc = 0) and 
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A"^'^^N{k = tt). For any given energy of the microcanonical 
ensemble, these quantities are computed as 



A"""n(fc = 0) 



J2a Waa{k = 0) 



=(fc = o)| 



iV, 



states "'mic 



f^mic(k — 0) 



and 



A'^'^Nik = tt) 



N states N„uc[k = tt) 



(7) 



(8) 



where naa{k — 0) and Naa{k ~ tt) are the eigenstate expec- 
tation values of n{k = 0) and N{k = tt), respectively, and 
nmic{k = 0) and Nmic{k = tt) are the microcanonical expec- 
tation values at any given energy E. The sum over a contains 
all states with energies in the window [E — AE, E + AE], 
and Nstates is the number of states in the sum. As discussed 
before AE = 0.1. 




FIG. 12: Average relative deviation of eigenstate expectation values 
with respect to the microcanonical prediction as a function of the 
effective temperature of the microcanonical ensemble. Results are 
presented for: (a) n{k = 0) and (b) N{k = tt). These systems have 
21 lattice sites, 7 bosons, and t — V = 1.0. An effective temperature 
T = 2 in these plots corresponds to an energy Emic ~ —5.25 for 
t' = V' ^ 0.03 and = -6.05 for t' = V ^ 0.24 and, T = 

10 corresponds to an energy Emic ~ —0.75 for t' — V' — 0.03 and 
Emic = —0.76 for t' = V' = 0.24 (see the corresponding parts of 
the spectrum in Fig.lllt. 

In order to compare results for different Hamiltonian pa- 
rameters, we have plotted A""=n(fc = 0) and A""=iV(fc ^ tt) 
as a function of the effective temperature T, corresponding 
to an energy E of the microcanonical ensemble. Having the 
energy of the microcanonical ensemble E, T can be com- 
puted from E = Z-^Tr^^Hcxp{-H/kBT)Y where Z = 

Tr|exp(-H/fcBT)}. Notice that here this effective temper- 
ature is only used as an auxiliary tool for assessing how far 
from the ground state these systems are, and to have a unique 



energy scale window for comparing different observables in- 
dependently of the Hamiltonian parameters t' and V', which 
change the ground state energy and level spacing. 

Results for A"*'=n(fc = 0) and A'^"'N{k = tt) are shown 
in Fig.[T2l'a) andfTSlb). respectively, for a system with L ~ 21 
and iVb — 7. The results depicted in Fig. [12] are qualitatively 
similar to those in Fig. 5 (main text). Below T ~ 1.5 fluctua- 
tions are in general large and nonmonotonic. For T > 1.5, one 
can see that A""'^ for n{k = 0) and N{k = tt) continuously 
increase, at any given temperature, as one approaches the in- 
tegrable point, i.e., there is no abrupt transition as t', V' 0. 

A direct comparison between some of the results presented 
in Fig. 5 (main text) and Fig.[T2]is depicted in Fig. [13] There, 
one can see that A""'^ in general reduces as one increases the 
system size [except for A"^^'^N{k = tt) when t' = V' = G\. 
The relative change is more pronounced away from integra- 
bility. These results do not allow us to discriminate between 
two possible scenarios for the breakdown of thermalization in 
quantum systems. A first scenario in which ETH may hold 
arbitrarily close to the integrable point as the system size in- 
creases, and another possible scenario in which in the ther- 
modynamic limit there may be a critical value of t', V' (the 
integrability breaking terms) below which ETH does not hold 
and above which it does. This question will need to be ad- 
dressed in future works. 




FIG. 13: Comparison between the average relative deviation of 
eigenstate expectation values with respect to the microcanonical pre- 
diction for two different system sizes. L = 21 and Nt = 7 (thick 
lines), and L = 24 and Nt — 8 (thin lines). Results are presented 
for: (a) n{k = 0) and (b) N{k = tt). In all cases t^V = 1.0. 
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